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ABSTRACT 

A new numerical technique is proposed for the electromagnetic characterization of 
the scattering by a three-dimensional cavity-backed aperture in an infinite ground plane. 
The technique combines the finite element and boundary integral methods to formulate 
a system of equations for the solution of the aperture fields and those inside the cavity. 
Specifically, the finite element method is employed to formulate the fields in the cavity 
region and the boundary integral approach is used in conjunction with the equivalence 
principle to represent the fields above the ground plane. Unlike traditional approaches, 
the proposed technique does not require knowledge of the cavity’s Green’s function and is, 
therefore, applicable to arbitrary shape depressions and material fillings. Furthermore, 
the proposed formulation leads to a system having a partly full and partly sparse as well 
as symmetric and banded matrix which can be solved efficiently using special algorithms. 



I. INTRODUCTION 


Recently, a new technique which combines the finite element and boundary integral 
formulations to yield a system for solution via the conjugate gradient method (CGM) 
and the fast Fourier transform (FFT) was proposed for electromagnetic scattering com- 
putations. In particular, the method was formulated for the characterization of filled 
slots and grooves in a thick ground plane with transverse electric (TE) and transverse 
magnetic (TM) incidence [1], [2]. A similar technique was also proposed in [3] for com- 
puting the aperture admittance matrix of the same geometry, but in that procedure a 
frontal solution algorithm was employed to reduce the memory demand. In this paper we 
describe an extension of the hybrid technique proposed in [1], [2] for the electromagnetic 
characterization of three-dimensional (3-D) cavity-backed apertures in a ground plane. 

The problem of scattering by 3-D cavity-backed apertures in a ground plane has 
been considered in the past via the mode matching [4] and moment method/modal [5] 
approaches. A unique feature of both of these techniques is a required knowledge of 
the cavity’s Green’s function. They are, thus, mostly restricted to rectangular cavities, 
but even then, additional difficulties may arise. For example, the mode matching tech- 
nique yields an infinite matrix that must be truncated and in the case of the moment 
method/moda! approach the admittance elements involve mode sums that are usually 
slowly converging. Neither of these methodologies are obviously capable of treating large 
size apertures but in the case of deep cavities, high frequency techniques such as those 
proposed in [6] and [7] could be effectively employed. Nevertheless, when the cavity is 
narrow in one dimension or shallow and perhaps filled with inhomogeneous material, a 
numerical solution approach is the likely alternative. 
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A numerical approach which has demonstrated promise for treating large structures 
is the finite difference- time domain approach (FD-TD) [8]. In this case the finite dif- 
ference mesh must be terminated with an absorbing boundary condition. In contrast, 
the proposed approach avoids possible truncation errors and excess discretization by 
representing the fields external to the cavity with the appropriate boundary integrals. 
Specifically, the proposed approach employs the finite element method [9] to formulate 
the fields within the cavity whereas the fields external to the cavity are expressed via 
the radiation integrals over the aperture. The resulting equations are then solved by 
demanding continuity of the tangential fields across the aperture. By virtue of the finite 
element method, the technique is applicable to cavities of arbitrary shape, possibly filled 
with inhomogeneous or composite materials. In the following sections the proposed hy- 
brid formulation is discussed in some detail along with the pertinent discretization of 
the resulting integral equations. A number of results are then presented which validate 
the formulation and we conclude with a discussion on the method’s merits and possible 
improvements. 


II. FORMULATION 

Consider the 3-D cavity-backed aperture illustrated in Figure 1. We will denote 
the free space region above the cavity (z > 0) as region I and that inside the cavity 
(— c < z < 0) as region II occupying the volume V. We will further assume that the 
cavity is filled with an inhomogeneous material having a relative permittivity c r (R) and 
relative permeability /x r (R). 

In accordance with the equivalence principle [10], the fields in the two regions can be 
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decoupled by closing the aperture with a perfect conductor and introducing the equiva- 
lent magnetic current 

M = E x z (1) 

over the extent of the aperture, where E is the electric field at the aperture (z = 0). The 
field in region I is then due to the radiation caused by the equivalent current M residing 
on the ground plane and possibly by other impressed sources (J 1 , M 1 ). Accordingly, by 
invoking image theory we have 

H J (R) = H‘(R) + H r (R) - 2jk 0 Y 0 JJ G 0 (R,R') • M(R ')dS' (2) 

where H* denotes the incident field due to (J*, M*) and H r is that reflected by the 
ground plane without the aperture. Also, ko = 2tt/A is the free space wavenumber, 
Zo = 1/Yq is the free space intrinsic impedance, S denotes the planar surface area of the 
aperture and Go is the free-space dyadic Green’s function given by 

S 0 (R,R') = (f + -|vv) G 0 {R,R') (3) 

with 

e -i*o|R-R'l 

I = xx + yy + zz and G 0 ( R, R') = 4?r | R _ R /| ■ 

Enforcing continuity of the tangential electric fields across the aperture, we find that 
the field in region II can be represented by the radiation of the equivalent magnetic 
current — M. The fields in the two regions are then coupled by enforcing continuity of 
the tangential magnetic fields across the aperture. This gives 

zxH J (M,J i I M i ) = zxH // (-M) z = 0 (4) 
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where H 1 and H 7/ denote the magnetic fields in regions I and II, respectively. 

Traditionally, M is solved from (4) by substituting for H 1 as given in (2) and ex- 
pressing as a function/integral of M. The resulting integral equation(s) are then 
discretized to obtain a system of equations for solution via direct inversion or LU de- 
composition. However, to explicitly express H n in terms of M implies a knowledge of 
the cavity’s Green’s function. For rectangular cavities filled with homogeneous material, 
this is usually found in modal form which is in general computationally inefficient [5]. 
Furthermore, in the case of arbitrarily shaped and/or inhomogeneously filled cavities 
there is no available closed form of the associated Green’s function. As a result, so far 
numerical solutions have only been considered for cavities that are rectangular and filled 
with homogeneous or uniformly layered material [4], [5]. 

To overcome the difficulty associated with the availability of the cavity Green’s func- 
tion, we employ the finite element method to formulate the fields in the cavity region 
(region II). Specifically, the cavity fields are demanded to satisfy the variational equation 

SF = 0 (5) 

where the functional F is given by [11] 

f, ( E ") ■ l#v[i( vxE ")'( v><E ") + ^( v,£ ' E "r 

-kle r E n • E 77 ] dV + jk 0 Z 0 if (e 77 x H 77 ) • MS (6) 

Scav 

if the variation is taken with respect to the electric field or by 

f ( H ") - 5j£[£( VxH '0*( VxH '0 + ^( v ** H ") J 

-klur H 77 • H 77 ] dV + jkoYo if (e 77 x H 77 ) • MS (7) 

5caw 
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if the variation is taken with respect to the magnetic field. In these, V denotes the 
volume occupied by region II, S cav corresponds to the surface that encloses V and n 
denotes the unit vector normal to the surface, pointing away from the cavity. Since 
the divergence of E 77 is zero in the source free region, it would appear that the second 
term in (6) and (7) is superfluous and, in fact, it is not included in the corresponding 
expressions found in [11]. However, a solution of (5) with 5 = 0 does not guarantee 
that the resulting fields will be maxwellian unless the divergence conditions are also 
satisfied. In the case of two-dimensional solutions they are satisfied a priori but not so 
for three-dimensional implementations. Here we propose that the divergence conditions 
be satisfied approximately and implicitly [12] by minimizing (6) or (7) with s £ 0. The 
parameter s is referred to as the penalty factor and its choice will be discussed later. 

To solve (6) or (7), it is necessary that the fields be known over the surface specified 
by S cav’ Obviously, the boundary conditions to be imposed on the conducting boundaries 
of the cavity are 

n x E /J = 0 (8) 

and ’ 

n • H /J = 0. ( 9 ) 

Substituting these into (6) and (7) eliminates the portion of the surface integral over 
the conducting boundary of the cavity (that is, there is no power flow through the 
metallic portion of S cav ). It remains to specify the boundary condition over the cavity’s 
aperture and this is given by (1). With these boundary conditions we may proceed with 
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a discretization of (5) which in conjunction with the discretization of (2) yields a system 
for a solution of M. Here the continuity condition (4) is enforced upon the discretized 
fields and this is the standard procedure employed in general finite element-boundary 
integral formulation [13]. However, if (6) is used the enforcement of (4) can be executed 
at the analytical stage by substituting (2) into (6). Doing so, we obtain the functional 


F = 


*// 


dV 


l III [i ( v * E ") • ( III. * V X E ") + ^ ( v • £ ' E ") 2 - *^ E " * E 

-2k 2 0 fj M(R) • ^jj s 5 0 (R, R') • M(R')dS'] dS 

-jk 0 Z 0 JJ s M(R) t [h’(R) + H r (R)] dS (10) 

which can be discretized via the finite element method for a solution of M. Once M is 
found, the far zone scattered field can be easily computed from 

H*(R) = - jk 0 Y 0 e -^- f( {66 + 4>4>) • M{x\y')e^^' sia ^dx'dy' (11) 

Lit XL JJS 

where ( R,9,4 > ) are the usual spherical coordinates of the observation point. The radar 
cross section (RCS) of the aperture is then given by 


“ llm 4,r * 2 i¥(ljr 


K— ► oo 


( 12 ) 


III. FINITE ELEMENT DISCRETIZATION 

For a numerical implementation of (5) we must first discretize the functionals by 
subdividing V and 5 into smaller volume and surface elements, respectively. Considering 
the electric field formulation, it is convenient to rewrite the functional F as 


F = Fy + Fs + Fe 


(13) 
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where Fy is the volume integral 


Fv =\\\i £ 


( 


dEz _ 9Ey\ 2 f dE x _ dEA 2 (dEy _ 8E, 

dy dz ) \ dz dx ) \ dx dy 


+■ 


»rt l T 


die^r) + die^) + 2 _ k 2 er ( E 2 + E 2 + E*) | dv ( 14 ) 


obtained by expanding the appropriate integrand in (10) and Fs denotes the surface 
integral 


F s = ~ 2k2 ° Us M ^ R) # [Us G °( R ’ R# ) • M ( R, )^ / 


dS. 


(15) 


The remaining portion of F is associated with the source field and is given by 


Fg = -2jkoZo 



M(R) • H’(R)d5 


(16) 


upon using the relation z x H r (R) = zx H*(R). For simplicity, in (14) we have omitted 
the superscript II and this practice will be continued in the remaining portion of the 
paper. 

To discretize (14) we subdivide the volume V into M v small volume elements such 
as tetrahedra, triangular prisms, or rectangular bricks. Within the eth element having 
n nodes, the field components are expressed as 

E P = Y1 Nj( x , y, z)<f>lj p = x,y,z (17) 

j'=i 

where N* are the known expansion or shape functions (see Appendix) chosen so that 4> pj 
(p — x , t/, z) represents the unknown field at the j th node of the eth element. Substituting 
(17) into (14) yields the portion of Fy attributed to the eth element. The complete 
expression for Fy is then obtained by summing/assembling the contributions from all 
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elements. This yields a functional in terms of the unknown node field components which 


must be found to satisfy (5). In accordance with the Rayleigh-Ritz procedure this is 
equivalent to setting the derivatives of F with respect to the node fields 4> p j ( p = x,y,z) 
equal to zero. Differentiating the portion of Fy attributed to the eth element with 


respect to the node field <j> e ri we obtain 


A fff 1 f \± djeJ^dj^NJ) dNf 

JJJv* Hr \ fr dx dx dy 


f dNj_ dNf dNf 
dy ^ dz dz 


s d(e r Nf) djc r Nf) dNf dNf 


V V ~ ' 

+ ej dx dz dz dx\* z ’\ 


dy dx 9y3 


Similarly, by differentiating with respect to the other node field components we have 


dFf_ _ y, fff 1 f [ a d(e r Nf) d(e r Nf) dNf dNf 

d<t> e yi JJJv* Mr \ c? dy dx dx dy 


\dNldNl . d( t ,N!)d(<rN') ar^Ml _ 1 

dx dx el dy dy dz dz 0 r r ' j yj 

Is dierNntKcrN*) dNldJl 11 
+ el dy dz ~ dz dy 


A fff 1 f [ d(e r Nn d(e r Nf) dNf 

JJJv‘ fi T t e.1 dz dx dx 


? dNf 1 

«* 


s d(e r Nf) d{e T Nf) dNf dNf ] dNf dNf 

el dz dy dy dz y: ^ dx dx 


We observe that if e r and fi r are assumed constant within the eth element, all integrals 
in (18)-(20) can be evaluated analytically. Otherwise, a numerical integration may be 


be required for their evaluation. 
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Let us now consider the discretization of the surface integral in (15). A difficulty 
in the evaluation of this integral is the usual singularity associated with the derivatives 
of the free space Green’s function. This, however, can be avoided by transfering the 
derivatives to the current. To do so, we invoke a common vector identity and the 
divergence theorem, leading to 


F s = - 2 Jf s M ( K )*{ k oJJ s M ( K ') G °( R ’ K ') dS ' 

+ V IIs W * M ( R, ) G °( R ’ R ') d5 ' 


dS. 


( 21 ) 


Through the same process, (21) can be further rewritten as 


F s = -2k\ JJ s M(R) • yj M(R/)G 0 (R, R ')dS' dS 
+2 V . M(R) yj s V' • M(R')G ? 0 (R, R !)dS' 


dS 


( 22 ) 


and by invoking (1) we obtain 


F S = -2/cq Jj s E x yj^ E x G 0 dS'j dS + 2k 2 0 jj^ E y E y G Q dS' 


dS 


(23) 


which can be discretized by subdividing S into M 3 smaller surface elements. In parallel 
with the volume discretization, the field components in the eth surface element can be 
expressed as 


n. 


e p = N !( x ^y)^ e P j p = x ’V 

i= 1 


(24) 


where n , denotes the number of nodes associated with the area element, N* are the 
same expansion functions as those in (17) with z = 0 and <f>*j (p = x,y ) represent the 
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node fields. Substituting (24) into (23) and replacing S in the first pair of integrals with 
S e , the area of the eth surface element, gives the portion of Fs attributed to the eth 
element. As noted earlier, to enforce the stationarity condition we need the derivatives 
of Fs with respect to the node fields. For the eth element we have 



and dF£/d<j> e zi = 0 since Fs is not a function of E z . We note that in deriving (25) and (26) 
the differentiation was performed only with respect to the node fields outside the square 
brackets in (23) while those introduced by substituting (2) into (6) remained uneffected. 
Further, we remark that the evaluation of the surface integrals in (25) and (26), although 
not trivial, can be done through analytical and numerical means as discussed in the 
Appendix. 

It remains to discretize (16) which involves the excitation fields. By replacing M 
with E in accordance with (1) we obtain 

Fe = 2jk 0 Zo Jj ( E X H * - E y H' x ) dS. (27) 

This can again be discretized by introducing the expansion (24) and by doing so we 
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obtain (for the eth element only) 


m 

m 


= 2jkoZ 0 Jf NfH*dS 


(28) 


and 


dF% 


= —2jkoZo ff N‘ H' x dS 


(29) 


Given the partial derivatives of all integral functions comprising the functional F we 
can now proceed with the construction of the final system of equations by imposing the 
stationarity condition (5). This implies that 


dF ^ dFt dF% ^ dF\ 


d<j>. 

dF 


- V' ZLV. + Y' , Y' _ o 

ZxW* 

MS , dFl , ^dF% 


dd> e • d<t> e • d<t> e • 

c— i e=i m/* e=i 


(30) 


dF _ %kdFl_ 

to* ~ htou 


i= 1,2, 3, ...,7V 


leading to a matrix system for the solution of the node fields. In (30) N denotes the 
total number of nodes, <f>pi ( p = z,y, z) are the node fields labeled with global indices 
and as before <p pi (p = x,y,z) are the node fields associated with the eth volume or 
surface element. Both <f> p i and <p pi refer to the field at the same node and thus the eth 
term of the summations have non-zero values only if the global node i belong to the 
eth element. The system implied by (30) must, of course, be solved after imposing the 
boundary condition (8) which permits us to zero out those field components that belong 
to nodes on metallic boundaries and are also tangential to that boundary. This reduces 
substantially the number of unknowns in the system which can then be solved via direct 
inversion, LU decomposition, or iteration. However, since the system matrix is partly full 
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and partly sparse as well as symmetric and banded (if the nodes are properly numbered), 
it can be more efficiently solved by those algorithms which exploit these properties [14]. 
Various partition techniques such as the those discussed in [15] can also be employed 
to enhance the efficiency of the solution. Further, the matrix system is amenable to a 
conjugate gradient-fast Fourier transform solution, thus, reducing the memory demand 
to 0(N) as in [1], [2], 

Let us now briefly address the formulation resulting from a finite element discretiza- 
tion of (7). In this case we cannot invoke (4) to incorporate (2) into (7) and obtain an 
expression that is completely in terms of the unknown nodal magnetic fields. As a re- 
sult, the continuity condition (4) must be explicitly imposed leading to a second system 
of equations to be coupled with that from the discretization of (7). The penalty with 
this approach is an increase in the number of unknowns because the surface magnetic 
currents or electric fields must be added to the nodal magnetic fields resulting from the 
discretization of region II. Accordingly, the system resulting from (7) is of the form 

[Km + [J?]{V>s) = o (31) 

where { <f > } is a column vector representing the nodal values of the discretized magnetic 
field in the cavity and over the aperture, and {^s} is a similar column for the discretized 
aperture electric fields. The matrix [ K ] is square and symmetric whereas [B] is rectangu- 
lar. Further, both matrices are very sparse and banded provided the nodes are properly 
numbered. A corresponding discretization of (2) via Galerkin’s methods leads to 

[B] T W + [P]{fe} = {4} (32) 

where the superscript T denote the transpose and [P] is a full, square and symmetric 
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matrix. Further, { b } is the excitation column vector and is a function of the incident field. 
Combining (31) and (32) yields a symmetric system that can be solved after imposing 
the boundary condition (9). 

IV. NUMERICAL RESULTS 

We present in this section some results for reference purposes and for validating the 
proposed formulation. In all cases, the excitation is a plane wave given by 

E’(R) = [(a • &)& + (a • W>' jk ’ ,R (33) 

where a = 0 X cos a + 4> x sin a is the polarization vector, k* is the propagation vector given 
by 

k* — — fc 0 ( sin 0 l cos <f> l x + sin 0 X sin <f>' y + cos O'i) (34) 

and O' and <j> % are the usual unit vectors in the spherical system and are associated with 
the angles O' and <£\ Further, in all computations the penalty factor s was set to unity 
and rectangular bricks were used for the discretization of region II. 

For reference, Figure 2 shows the y-component of the electric field along the center 
of an aperture formed by a narrow crack of length a — 0.5A, 1.0A, 1.5A and 2.0A. We 
remark that the behavior of the aperture field associated with the cavity-backed structure 
is distinctly different from that of a transmitting aperture as given in [16]. 

Figure 3 displays the co-polarized and cross- polarized backscatter RCS of a deep 
empty cavity as compared with data obtained via the moment method/modal approach 
[5], Overall the agreement between the results from the two methods is seen to be 
excellent except at the point {<j> = 90°) for a = 90°. We also note that the data in 
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this figure axe in agreement with those obtained using the mode matching procedure [4]. 
For the same aperture size, Figure 4 displays the cavity RCS as a function of depth. 
The resonant behavior of the cavity is rather distinct and characteristic of the cavity 
dimensions. 

The backscatter RCS of a material filled cavity is displayed in Figure 5. The re- 
sults based on this formulation and the moment method/modal approach are again in 
excellent agreement. Finally, Figure 6 refers to a cavity filled with multilayer material. 
The presented formulation is, of course, applicable to cavities filled with inhomogeneous 
material whereas traditional approaches are not. Further, arbitrary shape cavities can 
be treated with the same ease by employing an appropriate mesh generation algorithm. 

V. DISCUSSION 

There are a few issues which must be addressed in connection with the proposed 
formulation and the solution of the resulting system. In particular, below we discuss the 
role and choice of the penalty factor and the difficulty associated with the representation 
of the fields near corners and edges. Also, a few remarks are included with respect to 
the efficiency of this solution versus those based on more traditional approaches. 

A. Penalty Term 

As we stated above, the role of the penalty term is to implicitly satisfy the diver- 
gence conditions on the electric and magnetic fields. The term was first introduced in the 
stationary functional to remove the spurious modes that often appear in finite element 
solutions of eigenvalue problems such as those pertaining to waveguides and closed cavi- 
ties [12]. For those problems the role of the penalty term was very crucial in determining 
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the proper modes, propagation constant or resonant frequencies. For this application we 
observed that in most cases the actual value of the penalty factor ( s from 0 to 2) has 
a small effect on the RCS pattern. Nevetheless, it is desirable to avoid the ambiguity 
associated with the penalty term and this can be accomplished by using divergenceless 
basis functions [17]. 

5. Corner/Edge Conditions 

It is known [18] that some field components are singular near sharp perfectly conduct- 
ing edges as is the case with the four edges forming the aperture in Figure L The basis 
functions used in this paper or other similar conventional basis functions are not capable 
of accurately representing those field singularities. This, of course, leads to errors in the 
calculated field distribution which can be avoided by modifying the field expansion to 
include basis functions that simulate its expected singular behavior as determined from 
a static solution. A successful example of this is given in [19] and the approach could 
be extended to three-dimensions. We note that for the computations in this paper we 
did not observe a noticeable effect on the RCS pattern due to errors associated with the 
edge condition. 

C. Computational Efficiency 

The proposed formulation was considered because of its potential to treat cavities of 
arbitrary shape and material fillings. However, it was also found more computationally 
efficient than tradiational formulations [4], [5] without even making use of the symmetry 
and sparseness properties of the resulting matrix. In particular, by using LU decom- 
position for the solution of the matrix system, the data in Figure 3 were generated in 
less than one minute on an Apollo DN10000. In contrast, the corresponding time when 
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using the mode matching technique was five minutes on a Cray XMP-48 [4] and that 
for the moment method/modal approach was two minutes on an Apollo DN10000. The 
CPU time differences between the moment method/modal approach and the proposed 
formulation are even more apart for larger apertures. For example, the data in Figure 
5 were obtained in four minutes on an Apollo DN10000 by this method but 65 minutes 
were required with the moment method/modal approach on the same machine. A major 
reason for the large difference in CPU time is due to the slowly converging mode sums 
that must be evaluated in the process of generating the matrix elements. 

VI. SUMMARY 

A new technique was proposed for a numerical characterization of the scattering 
by a 3-D cavity-backed aperture in a ground plane. The proposed technique combines 
the finite element and boundary integral methods to generate a system of equations 
for the aperture fields in conjunction with the continuity condition at the aperture. 
In principle, the technique is capable of treating arbitrarily shaped cavities filled with 
inhomogeneous materials. However, because of the lack of available reference data, the 
proposed solution was only validated for rectangular cavities. Some important issues 
relating to the implementation and effciency of the technique were also discussed. 

APPENDIX 

A suitable element for the discretization of a rectangular cavity is the rectangular 
brick illustrated in Figure 7. For this element, the expansion functions in (17) are given 
by 

{£ - xjjd - yW_ - *0 x'W - yQ(c' - d) _ * W - jO 

1 a'b'cf 2 a'Vd 3 a'Vc' 
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( 35 ) 


NZ 

Nj 


(g[ - *'¥0' - *0 iye = ( fl/ - - gV 

a'6'c' 5 a'ftV 

£V£ jy e («' - *'¥*' 

aW 8 a'Vc' 


x'jbf - yV 

a'b'c' 


and those in (24) are obtained by setting z f equal to zero. Here (x\y f ) denote the local 
coordinates associated with the e element and through a linear transformation they can 
be replaced by global coordinates. Substituting these into (18)-(20), (28) and (29), the 
resulting integrals can be evaluated analytically on the assumption of constant e r and ji T 
within each element. The evaluation of the integrals in (25) and (26) is, however, more 
involved because of the kernel’s singularity. To illustrate how these integrals could be 
evaluated let us consider one that has the generic form 


Pi? = jj st N!(z,y) [Jl N?(x\y')GoCR,R')dS’ 


dS 


(36) 


If the eth and e'th elements are not adjacent, then mid-point integration could be used 
to evaluate p\f as 


Pit = 


(37) 


where the subscript c denotes the point at the center of the element. Otherwise, the 
eth and e'th elements could be subdivided into 3x3 small rectangles. Replacing the 
expansion functions Nf and Nf within each of nine rectangles with their mid-point 
values yields 


Pij — ^2 ^2 Vmn)Nj ( ^pq^V^Qmnpq (38) 

m,n= 1 p 7 q—l 

where 


Qmnpq — 



<7o(R,R')<H" 


dS. 


(39) 
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The integral g m npq can again be evaluated via mid-point integration provided (x,y) and 
( x f ,y f ) do not belong to the same rectangle, or share the same side. Otherwise, we 
rewrite g mnp q as 


9mnvq ~ 4x IL list 


,-jko |R-R'| 


|R - R'| |R — R'| 


dS' 


dS 


+ 


7" // ( II ip~ 1 W 7\ dS '] dS - 

4x JJs \JJs‘' q |R - R | J 


(40) 


The first integral of these has a non-singular integrand and can therefore be evaluated 
numerically using a 2 x 2 Gaussian integration. The second integral has, of course, a 
singular integrand and must be evaluated analytically as described in [5]. 
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FIGURE CAPTIONS 


Fig. 1 Geometry of a cavity-backed aperture in a ground plane. 

Fig. 2 Aperture field distribution ( E y ) at the center of a narrow crack, b = 0.05A, 
c = 0.25A, e r = 1.0, Hr = 1.0, normal incidence, (a) a = 0.5A. (b) a = 1.0A. (c) 
a = 1.5A. (d) a = 2.0A. 

Fig. 3 Backscatter RCS patterns for an empty cavity versus incidence angle, a = 0.7A, 
b = 0.1 A, c = 1.73A, 6 = 40°. Solid and dashed lines represent the solution of this 
method; circles and squares represent the moment method/modal solution [5]. (a) 
a = 90° (E = <££*). (b) a = 0° (E = BE$). 

Fig. 4 Backscatter RCS for an empty cavity as a function of cavity depth, a = OJA, 
b = 0.1A, 6 = 40°. 

Fig. 5 Backscatter RCS patterns for a material-filled cavity versus incidence angle. 
a = 1.0A, b = 0.25A, c = 0.25A, e r = 7.0 — j0.5, n T = 1.8 - jO.l. Solid and 
dashed lines represent the solution of this method; circles and squares represent 
the moment method/modal solution [5]. (a) <f>= 90°. (b) <j> = 0°. 

Fig. 6 Backscatter RCS patterns for a multilayer material-filled cavity versus incidence 
angle, a — 1.0A, b = 0.25A, c = 0.25A. Top layer: e r = 7.0 — j’0.5, \x T — 1.8 — jO.l, t 
(thickness)= 0.0625A; middle layer: e r = 3.0- j0.05, = 2.0, t = 0.125A; bottom 

layer: e r = 5.0 — j0.03, fi T = 1.0, t = 0.0625A. Solid and dashed lines represent the 
solution of this method; circles and squares represent the two-dimensional solution 
[1], [2]. (a) 4> = 90°. (b) 4> = 0°. 


Fig. 7 Geometry of a rectangular brick. 
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